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Abstract 

We assess the risk of an Earth impact for asteroid (99942) Apophis by means 
of a statistical analysis accounting for the uncertainty of both the orbital so- 
lution and the Yarkovsky effect. We select those observations with either 
rigorous uncertainty information provided by the observer or a high estab- 
lished accuracy. For the Yarkovsky effect we perform a Monte Carlo simula- 
tion that fully accounts for the uncertainty in the physical characterization, 
especially for the unknown spin orientation. By mapping the uncertainty in- 
formation onto the 2029 fe-plane and identifying the keyholes corresponding 
to subsequent impacts we assess the impact risk for future encounters. In 
particular, we find an impact probability greater than 10~ 6 for an impact in 
2068. We analyze the stability of the impact probability with respect to the 
assumptions on Apophis' physical characterization and consider the possible 
effect of the early 2013 radar apparition. 
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1. Introduction 

Asteroid (99942) Apophis was discovered by R. A. Tucker, D. J. Tholen, 
and F. Bernardi at Kitt Peak, Arizona on 2004 June 19 (Minor Planet Sup- 
plement 109613). Because of its motion, the asteroid was immediately recog- 
nized as an object of interest and therefore was observed again on June 20. 
In the following days, issues related to telescope scheduling, bad weather, 
and lunar interference made it impossible to collect new observations. By 
the time the Moon was out of the way, the solar elongation had become too 
small, thus making the object even more difficult to observe. Thus, the aster- 
oid was lost until December 2 004, when it was redis covered by G. J. Garradd 



at Siding Springs, Australia (IGilmore et al 



20041 ). Thanks to the Decem- 



ber 2004 observations, Apophis was recognized as a potentially hazardous 
asteroid, with a peak impact probability (IP) of 2.7% in April 2029 reported 
by Sentrjtl and NEODySB Also, the new observations revealed astrometric 
errors in the discovery observations, for which D. J. Tholen remeasured the 
positions and adjusted the time (Minor Planet Circulars 53585 and 54280). 



On 2004 December 27 t 



vat ions ( IGleason et al 



re Sp acewatch survey reported pre-discovery obser- 



20041 ) from March 2004 that ruled out any impact 



possibility for 2029. 

In January 2005 Apophis was observed by radar from the Arecibo ob- 



1 http://neo.jpl. nasa.gov/risk 
2 http: / / newton.dm.unipi.it / neodys 
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servatory, obtaining three Doppler and two delay measurements. In August 
2005 and May 2006 two more Doppler observations were obtained. These 
observations significantly reduced Apophis orbital uncertain ty and led to a 



6 Ea rth radii nominal distance from the geocenter in 2029 ( jGiorgini et al 



20081 ). Though the 2029 impact was ruled out, other potential impacts were 



revealed in the following decades. Moreover, the exceptionally small close 
approach distance turns a well determined pre-2029 orbit to a poorly known 
post-2029 orbit for which scattering effects are dominant. Therefore, even 
sm all pertu r bation s prior to 2029 play an importan t role. 



Chesleyl (120061 ) shows that the Yarkovsky effect (IBottke et al.l . 120061 ) sub- 



stantially affects post-2029 predictions, and therefore h a s to 
account for Apophis impact predictions. 



3e taken into 



Bancelin et al. 



( 120121 ) also suggest 



that the Yarkovsky effect is relevant, though they do not account for the 
corresponding perturbation for their risk assessment. However, the avail- 
able physical characterization of Apophis is not sufficient to estimate the 
Yarkovsky effect as was done in the cas e of asteroid (101955) 1999 RQ36 



( IMilani et al. 



2009 



Chesley et al.l 



20121 ). In particular the lack of infor- 



mation for spin orientation does not allow one to determine whether the 
Apophis semimajor axis is drifting inwards or outwards. In some cases, 
the orbital drift associated with th e Yarkovsky effect can b e measured from 



the orbital fit to the observations ( Farnocchia et al. 



20131)- Unluckily, this 



is not yet possible for Apophis, and upcoming radar observations are un- 
likely to help in constraining the Yarkovsky effect. Another nongravitational 



perturbation that coulc 



je im portant is solar radiation pressure. However, 



Zizka and Vokrouhlickyl ( 1201 ll ) show that solar radiation pressure has a much 
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smaller effect on the Apophis trajectory than does the Yarkovsky effect. 



2. Orbit determination 



2.1. Dynamical model 

The scattering close approach in 2029 calls for a dynamical model that 
is as accurate as possible. We used an N-body model that includes the 
Newtonian attraction of the Sun, the planets, the Moon, and Plut o. These 



acce 



erations are based on JPL's DE424 planetary ephemerides ( jFolkner 



201 lh . We further added the attraction of 25 perturbing asteroids, which are 
listed in Table ffl 

As was already pointed out by 



Chesley et al. 



(12012( 1 and 



Farnocchia et al 



( 120131 ). the relativistic terms of the planets are important when dealing 



with high precision orbit modeling. Th erefore, we used the Einstein-Infeld- 
Hoffman approximation ( iMoyerl . 120031 ). which accounts for planetary rela- 
tivistic effects. 



2.2. Astrometry 

We aimed at u sing only high qualit y astrometry. Therefore, we selected 



observations from 



Tholen et al. 



(120121 ). for which rigo rous uncertai n ty in 



formation is provided by the observer. In particular, 
quantify the error due to different sources: 

• astrometric fit error Aa a and A5 a ; 



Tholen et al. 



torn 



centroiding error Aa x and A5 



tracking error Aat and A5t- 
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These observations were obtained by Mauna Kea an d Kitt Peak observato - 



ries and were reduced using the 2MASS star catalog (ISkrutskie et al. 



20061 ) 



Though 2MASS has been recently u sed as a reference to remove system- 



atic e rrors present in other catalogs (jChesley et al 



2010) 



Zacharias et al 



(120031 ) discuss the presence of biases in the 2MASS cata 



tematic error s in de clination are reported by 



Tholen et al. 



Tholen et al. 



og, an d 0.02" sys- 
(120 lj ) . Moreover, 



( 120131 ) show how the absence of proper motion in the 2MASS 



catalog is responsibl e for regional biases detected in 2011 Pan-STARRS PS1 



asteroid astrometry (IMilani et al. 



20121 ). The bias in declination due to the 



absence of proper motion is fairly close to a mean value 0.05". On the other 
hand the right ascension bias is variable, getting as large as 0.14", and shows 
a strong regional signature. Even so, 0.05" is a reasonable mean for the right 
ascension bias too. Therefore, to mitigate the absence of proper motion, we 
relaxed the observational weights by adding a time dependent uncertainty 
component related to proper motion biases according to 

yr - 2000 



a PM = 0.05" 



(1) 



2011 - 2000 
where yr is the year of the observation. 

This time-dependent component of the bias takes into account the uncor- 
rected proper motion between the epoch of the observation and the reference 
epoch of the 2MASS catalog (J2000). However, since the 2MASS catalog 
has been obtained from images taken over a period of about 3 years, and not 
instantaneously at the reference epoch, even the nominal J2000 positions in 
the catalog may be biased by an amount directly related to how long before 
the year 2000 they were actually observed by the survey. Therefore, we con- 
servatively assumed an additional <jb — 0.04" error component that accounts 



for the intrinsic bias of the catalog at the J2000 epoch. 



Tholen et al 



(120 12[ ) 



By rolling up all the sources of error, we gave the 
observations the following individual weights: 

w a = l/(Aa 2 a + Aal + Aa 2 t +Nal + Na 2 PM ) (2) 
w 5 = l/(A5 2 a + A5 2 x + A5 2 + Na 2 B + Na 2 PM ) (3) 

where N is the number of observations in the same night and accounts for the 
fact that catalog biases of observations obtained in the same region of the sky 
are expected to be highly correlated (since the same star s wer e used in the 



astrometric solution). The results in 



Milani et al 



(120121 ) and 



Tholen et al 



(120131 ) show that the angular scale of the correlation in the proper motion 
bias is larger than the typical sky-plane motion of Apophis in one night. 
Therefore, the temporal correlation of the biases likely extends to more than 
a single night of observations, and N should be increased to include all ob- 
servations within a certain angular distance of each other (up to a few tens 
of degrees, and location-dependent). A better statistical treatment would be 
to correct the observations by subtracting the mentioned biases. However, 
this is beyond the scope of the present paper and we plan to address this 
issue in the future. 

In our analysis we also include a 2MASS-based position from the March 
2004 precovery observations by the Spa cewatch survey, rem easured with the 



Tholen et al. 



(120121 ). from a stack of 



same techniques and tools described in 
five of the original exposures. Formal uncertainty values are also available 
for this position, and are treated in the same way. 

Apophis was observed by the Arecibo radar telescope in 2005 and 2006. 
This radar apparition produced 2 delay and 5 Doppler measurements for 
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which rigorous uncertainty information is available ( jGiorgini et all 120081 ). 

The astrometry included so far covers the time interval from 2004 March 
15 to 2008 January 9. However, Apophis has been further observed through 
late 2012. Therefore, we included observations from Magdalena Ridge (four 
observations on 2011 March 14) and Pan-STARRS PS1 (four observations 
on 2012 August 8 and four on 2012 December 29). These are the only two 
observatories with a residual RMS for numbered asteroid observations below 
0.2" according to AstDyqfi Since the Pan-STARRS PS1 observations are 
reduced using the 2MASS star catalog, we relaxed weights consistent with 
(ED and (ED: 



w a = 1/(0.30" 2 + No~ 2 B + Na 2 PM ) 
w 5 = l/(0.30 ,,2 + Na 2 B + Na 2 PM ) 



(4) 
(5) 



On the other hand, Magdalena Ridge observations are reduced using the 



USNO-B1.0 catalog (IMonet et al. 



We therefore debiased according to 



2003) , which a c count s for proper motion. 



Chesley et al. 



( 120101 ) and weighted these 



observations at 0.30". These twelve observations doubled the observed arc 
with quality comparable to the 



Tholen et al. 



(120121 ) astrometry. 



2.3. Orbital solution 

By using the observational data described in Sec. 12.21 we obtain the os- 
culating orbital elements in Table |2j We t hen mapped the orbita l elements 
and their uncertainties to the 2029 6-plane (jValsecchi et all 120031 ) . Figured] 
shows the positions on the 6-plane corresponding to the orbital solutions: T 



5 http://hamilton. dm. unipi.it/astdys2 
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is for only Tholen et al. observations, TR is for Tholen et al. plus radar 
observations, and TRO is for Tholen et al., radar, plus Pan-STARRS PS1 
and Magdalena Ridge observations. 

Table |3] reports 6-plane coordinates and uncertainties for different choices 
of the observational dataset. (Note that the 6-plane is normal to the asymp- 
totic velocity and the (£, () coordinates mark the intersection of the asymp- 



tote with respect to the geocenter. Therefore, b = \f£, 2 + C 2 does not cor- 
respond to the closest approach, which is at about 39000 km.) We also 
included the 6-plane coordinates for the orbital solution obtained by using 
all the available observations. In some cases there was an excess of observa- 
tions from a single observatory in a single night. In order to reach the level 
of a preferred contribution of 5 observations per night we relaxed the weights 
in these cases by a factor a/ N/5, where N is the number of observation in 
the night from the same observatory. 

It is worth noting that the T and TRO nominal solutions are signifi- 
cantly shifted with respect to the TR prediction: the distances in a-space 
are 1.85 and 1.46, respectively. This difference may point to the need for a 
better data treatment that consists in correcting observations for catalog and 
proper motion biases, as discussed in Sec. 12.21 It is also interesting to note 
that the solution obtained with the full dataset is at 0.60 a from the TRO 
prediction and the impro vement in the orbital uncertainty is about 20%. As 



Krolikowska et al. 



( 120091 ) show that an appropriate selection of the astromet- 
ric weighting scheme may be essential for impact predictions, we prefer to 
use TRO as nominal solution, as we have observations with a more reliable 
internal consistency and uncertainty estimate. 
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Figure 1: Nominal 2029 6-plane coordinates and 3-ct uncertainty ellipses for different sets 
of observations: T is for only Tholen et al. observations, TR is for Tholen et al. plus radar 
observations, and TRO is for Tholen et al., radar, plus Pan-STARRS PS1 and Magdalena 
Ridge observations. Note that the axes do not use the same scale. 

3. The Yarkovsky effect 

The Yarkovsky effect is a nongravitational perturbation that causes aster- 
oids to undergo a secular variation in semimajor axis, resulting in a quadratic 
in time runoff in a nomaly. 



As described in 



Farnocchia et al. 



( 120131 ). we modeled the Yarkovsky effect 



as a purely transverse acceleration a t = A 2 /r 2 , where r is the heliocentric 
distance and A 2 is a function of the asteroid's physical quantities. The cor- 
responding semimajor axis drift can be computed as 

da _ 2A 2 (l-e 2 ^ 
dt np 



, • (6) 
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where e is the eccentricity, n is the mean motion, and p is the semilatus 
rectum. 

As mentioned, A2 can be neither computed from the available physical 
characterization nor estimated from the orbital fit. In fact, with current 
astrometry we get an uncertainty <ja 2 = 70 x 10~ 15 au/d 2 while the likely 
range of A2 for an object of this size is ±25 x 10~ 15 au/d 2 . Therefore, to 
compute the distribution of the semimajor axis we performed a M onte Carlo 



simulation as was previously suggested by 



Chesley et al. 



(120091 ). Assumed 



physical parameters for the simulation were obtained as follows: 

• Recent observations from ESA's Herschel space observatorja give a 
diameter D = (325 ± 15) m (T. Mueller, personal communication). 
Thus, we used a normal distribution with mean 325 m and standard 
deviation 15 m (see top left panel of Fig. [2]). Herschel observations 
also provide a geometric albedo py = 0.23 ± 0.02. By assuming a 
slope parameter G = 0.1 5 ± 0.05, we c o nvert ed py to the Bond albedo 
A = (0.29 + 0.684G)p v (IBowell et all . Il989l ). The top right panel of 
Fig. [2] shows the distribution of A. 



Binzel et al. 



( 120091 ) report a grain density between 3.4 and 3.6 g/cm 3 
and a total porosity between 4% and 62%. This leads to a density 
p between 1.29 and 3.46 g/cm 3 . Therefore, we assumed a lognormal 
distribution with mean 2.20 g/cm 3 and variance 0.3 g 2 /cm 6 (see middle 



lef t panel of Fig 



of 



Binzel et al. 



. This distribution is consistent with the constraints 
( 120091 ). as the probability of being smaller than 1.14 



1 http://www.esa.int/Our_Activities/Space_Science/Herschel 
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g/cm 3 or larger than 3.24 g/cm 3 is 2.03% and 2.47%, respectively. 



Delbo et al. 



For Ap ophis the thermal inertia has not yet been measured. 
( 12007b! ) find a generic relationship between the thermal inertia and the 
diameter, according to which F = doD~^, with do = (300 ± 45) J m -2 



-0.5 



K 



D is in km, and ip = 0.36 ± 0.09. Accordingly, we used 



normal distributions for F and tp, and combined them with the already 
available diameter distribution. The middle right panel of Fig. [2] shows 
the distribution of the thermal inertia. 

Behrend's asteroids and comets rotation curves websitej reports a ro- 
tation period P rot = (30.4008 ± 0.0144) fjj, therefore we used the cor- 
responding normal distribution (bottom left panel of Fig. [2]). 

The obliquity of the Apophis rotational pole is unknown. Assuming a 
uniform distr ibution for the o bliqui ty is a poor approximation. In fact, 
according to lLa Spina et al.l (120041 ) retrograde and direct rotators are 
in a 2:1 ratio within the NEO population. As shown in the bottom right 



panel 



of Fig. [21 we used the obliquity distribution from 



Farnocchia et al 



( 120131 ). where the obliquity is inferred from a list of detections of the 



of 



Yarkovsky effect. 

For each poi nt in the Monte Carlo s ampling we computed the associated 

(120131 ). which is based on the linear model 



A ? according to 



Farnocchia et al. 



Vokrouhlicky et al. 



( 120001 ). and the orbital drift was then obtained by ()6j). 



5 http://obswww.unige.ch/^behrend/page_cou.html 
6 http:/ /obswww.unige.ch/^behrend/r099942a.png 
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Figure E] shows the distribution of A 2 and da/dt. The bimodality of the 
distributions corresponds to the extreme obliquity peaks. 

As previously done for the orbital solution uncertainty, we want to see 
how the Yarkovsky effect maps to the 2029 6-plane. To do so we assumed 
an arbitrary reference value (A2) re f = 10 x 10~ 15 au/d 2 , corresponding to 
(da/dt) re f = 4.59 x 10~ 4 au/Myr, for which we obtained A( re f = 138.18 
km and A£ re / = 0.45 km. This confirms that the Yarkovsky perturbation 
essentially affects only th e ( coordinate, which is strongly related to the time 
and then to the anomaly ( IValsecchi et al.1 120031 ). To compute the value of A£ 
for a generic point of the Monte Carlo sampling we scaled from the reference 
value, after testing the validity of the linear approximation. 



4. Keyholes and risk assessment 



Even though the 2029 close approach does not lead to an i mpact, we can 



Valsecchi et al 



use th e corresponding 6-plane analysis to predict future impacts. 
( 120031 ) show that the 6-plane coordinates leading to a resonant impact lie on 
predictable circles. Figure H] shows the Valsecchi circle on the 2029 6-plane 
for the 2036 reso nant return, c orrespo nding to a 7:6 resonance. The Line of 



Variation (LOV, 



Milani et al. 



2005bl ) is representative of the orbital uncer- 



tainty region. The intersecti on between the orbital uncertainty region and a 
Valsecchi circle is a keyhole ( jChodasl . Il999l ). leading to a future impact. 

To find the keyholes on the 2029 6-plane we propagated a densely sam- 
pled LOV through the year 2113, recording the circumstances of all close 
approaches to Earth and in particular the 6-plane coordinates, i.e., (£, £). 
On the 6-plane of a given post-2029 encounter we can interpolate among 
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nearby samples to identify the minimum possible encounter distance along 
the LOV. When this minimum distance was smaller than 1 Earth radius 
(-Re), we obtained the keyhole width by mapping back to the 2029 6-plane 
the chord corresponding to the intersection between the LOV and the Earth 
cross section. This technique allows us to identify keyholes in 2029 as small 
as a few centimeters in extent (Table H]). If we plot the minimum encounter 
distance for each LOV sample over the 100 year integration (and each inter- 
polated local minimum in the encounter distance) against the associated 2029 
6-plane coordinate C2029, then we obtain the plot shown in Fig. [5j The upper 
panel shows the wider view, which has been extended to the left from nom- 
inal to capture the complex of secondary potential impacts surrounding the 
2036 primary resonant return that is situated at about —1500 km on the ab- 
scissa. There are other similar collections of potential impactors in the plot, 
but the 2036 case is the only one in the depicted ( range where the primary 
return can lead to impact. As an example, the lower panel in Fig. [5] pro- 
vides a detailed view of the complex of potential impactors situated around 
—220 km on the plot. Here the central post-2029 return, which spawns the 
surrounding secondary returns, is in 2051, but no impact is possible in that 
year, although an approach closer than 2i? e is possible. 

Figure [6] shows the ( distribution on the 2029 6-plane by combining 1) the 
orbital uncertainty obtained from the orbital fit to the observations neglecting 
the Yarkovsky effect, and 2) the dispersion due to the Yarkovsky effect. 
The overall ( distribution is therefore obtained as the convolution of the 
starting probability density functions. It is interesting to see how the overall 
probability density function (PDF) is close to the Yarkovsky PDF. This is 
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due to the low dispersion of ( related to the orbital uncertainty PDF: if 
the latter were a Dirac delta function, the overall PDF would coincide with 
the Yarkovsky PDF. The figure also contains the keyholes, where the height 
of the corresponding bars is determined by the keyhole width. The impact 
probability can be computed as the product of the width of the keyholes and 
the PDF. Table H] reports the information on the keyholes along with the 
associated impact probability. 

The results contained in Table H] depend on the assumptions made on the 
physical modeling. For example, using a distribution for which a retrograde 
rotator is so much more likely than a direct one could be inaccurate. In 
general, we expect that the impact probabilities associated with keyholes 
close to the peaks of the PDF should be more stable, while those associated 
with keyholes in the tails of the distribution should be more sensitive to 
changes in the assumptions because of the exponential decay. To confirm 
this idea, Table reports the values of the impact probability of the 2036 
and 2068 (A£ = —212.14 km) keyholes, for different assumptions on Apophis' 
physical characterization. We can see that the order of magnitude of IP2068 
does not change, while IP2036 changes up to three orders of magnitude. 

The e xceptionally clo s e app roach of 2029 may change Apophis' spin ori- 



entation. 



Scheeres et al. 



(120051 ) show that Apophis could exit the 2029 en- 
counter in complex rotation. Therefore, the post-2029 dynamics is a function 
°f C2029 and post-2029 A 2 , which we denote with A' 2 . In particular, the 2029 
keyholes should be considered a two-dimensional entity in (C2029, -42)- Then, 
the IP of a given keyhole is 

/+00 
(h(( W 29)w) f 2 (A 2 )dA' 2 (7) 
-00 
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where f\ and f'2 are the PDF of C2029 and A' 2 , and w is the keyhole width. 
We checked that for different values of A' 2 the keyhole location changes only 
slightly (up to 0.5 km) and the width does not change significantly. Therefore, 
IP does not depend significantly on A' 2 . 

5. Early 2013 apparition 

As already mentioned, the observational dataset available for Apophis 
does not allow us to measure the Yarkovsky effect from the orbital fit. 
Apophis is scheduled to be observed from the Goldstone radar observatory in 
January 2013] and the Arecibo radar observatory in February 2013o There- 
fore, we wondered whether or not this radar apparition will constrain the 
Yarkovsky effect. To answer this question we considered three different sce- 
narios: 

(a) no Yarkovsky effect, i.e., A2 = au/d 2 ; 

(b) positive orbital drift with A 2 = 25 x 10 -15 au/d 2 , which corresponds 
to the right peak of the da/dt distribution of Fig. EJ 

(c) negative orbital drift with A2 = —25 x 10~ 15 au/d 2 , which corresponds 
to the left peak of the da/dt distribution of Fig. [3j 

Figure [7] shows the predicted values of delay and Doppler measurements 
for 2013 February 6, for the different scenarios of the Yarkovsky effect. The 
uncertainty ellipses overlap each other and their intersection contains all the 
nominal predictions. We conclude that the Yarkovsky effect will not be 



7 http://echo.jpl. nasa.gov/asteroids/Apophis/ Apophis_2013_planning.html 
8 http://www. naic.edu/ ^pradar/sched.shtml 
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detectable with the next radar apparition. To further quantify this, we sim- 
ulated synthetic delay measurements on 2013 February 5 and 6 (see Table [6]). 
As the relative accuracy of delay measurements is ~100 times better than the 
relative accuracy of Doppler measurements, two delay measurements provide 
a stronger orbital constraint than a single dela y/Doppler pair. For each sce- 



nario we tried to determine A<i as was done in iFarnocchia et al.l (120131 ) . We 
found <ta 2 — 48 x 10~ 15 au/d 2 with the simulated radar data, which is about 
two times larger than a reasonable value of \A 2 \. 



This result may look surprising, s ince in the cases o 



f lCheslev et all 120031 ) and 1999 RQ 36 flCheslev et all 120121 ) the availability 



(6489) Golevka 



of three radar apparitions allowed strong detections of the Yarkovsky effect. 
However, in these cases the time span covered by radar observations was 12 
yr (1991 to 2003 for Golevka, 1999 to 2011 for 1999 RQ 36 ) while for Apophis 
it is 8 yr. Moreover, the second radar apparitions of Golevka in 1995 and 
1999 RQ36 in 2005 were approximately in the middle of the other two, while 
the second apparition of Apophis was in 2006, just one year after the first 
apparition, and contains only one Doppler measurement, which is not as 
precise as delay measurements. Therefore, both the timing and precision of 
the data may explain the failure in detecting the Yarkovsky effect with three 
radar apparitions. As a matter of fact, the 2005, 2006, and 2013 datasets 
essentially contain the same information as two radar apparitions, since only 
2005 and 2013 contain delay measurements. 

At this point it is interesting to see what would be the effect of the 
January/February 2013 radar apparition on the risk assessment. As the 
Yarkovsky effect will not be directly detected, the main contribution will be 
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on the orbital uncertainty information. For example, let us consider scenario 
(c), i.e., we used the radar observations simulated according to the assumed 
A2 = —25 x 10~ 15 au/d 2 . We computed the corresponding orbital solution 
without accounting for the Yarkovsky effect, which is not measurable. The 
corresponding 2029 6-plane coordinates are (£,C) — (9485.10,47369.09) km, 
with uncertainties 5.30 km and 11.00 km, respectively. While the uncertainty 
in £ is almost the same, the uncertainty in ( decreased by a factor of 4 and 
then the corresponding PDF of ( due to the orbital uncertainty is even closer 
to a Dirac delta function. As a consequence, the Yarkovsky PDF and the 
overall PDF have almost the same curve. The changes in the IP are mostly 
due to a shift of about -3 00 km of the nomin al Yarkovsky-free prediction. 
All the virtual impactors ( IMilani et all |2005aj) of Table H] persist, and four 
of them now have an IP> 8.6 x 10~ 8 : IP 20 36 = 2.23 x 10~ 7 , IP 2 o6s(AC = 
-212.14W) = 6.95 x 10~ 7 , IP 206 8(AC = 47.68A;m) = 1.03 x 10~ 7 , and 
IP 2105 = 2.79 x 10- 7 . 

Though the Yarkovsky effect will not be constrained, radar and optical 
photometry can provide independent information to improve Apophis' phys- 
ical modeling. In particular, it may be possible to further constrain the 
diameter and the albedo, and to estimate the spin orientation, which would 
allow us to remove the bimodality in the Yarkovsky distribution. 

The next opportunity for radar observations is the September 2021 close 
approach. Meanwhile, ground based optical observations are possible and 
can help in determining the Yarkovsky effect before the September 2021 
close approach. Table [7] contains the uncertainty in the A 2 determination by 
assuming simulated radar observations in February 2013 and progressively 
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adding simulated ground based optical observations weighted at 0.15". We 
conclude that the Yarkovsky effect would be significantly constrained between 
late 2020 and early 2021. 

6. Conclusions 

In this paper we presented an updated and careful impact risk analysis for 
asteroid (99942) Apophis. The presence of the scattering close approach in 
2029 calls for accurate orbit estimation and dynamical modeling. Therefore, 
we selected astrometry for which the uncertainty information is rigorously 
computed by the observer and added the contributions of observatories with 
comparable performance. To model the Yarkovsky effect we performed a 
Monte Carlo simulation relying on the best physical model currently avail- 
able. 

We obtained the distribution on the 2029 6-plane by combining the or- 
bital uncertainty and the dispersion due to the Yarkovsky effect. We then 
found twenty keyholes corresponding to future impacts and computed the 
corresponding impact probabilities. In particular, we found a 2.27 x 10 -6 
impact probability for an impact in 2068. 

We analyzed the stability of the computed impact probabilities with re- 
spect to the assumed Apophis' physical characterization. Impacts corre- 
sponding to keyholes close to the nominal solution have a stable impact 
probability (the order of magnitude does not change). On the other hand, 
for keyholes in the tails of the distribution we found variations up to three 
orders of magnitude. 

We also showed that the February 2013 radar apparition will not allow 
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a direct detection of the Yarkovsky effect, and therefore the post-2029 pos- 
sibility of impact will not be ruled out. The main contribution will be a 
reduction of the orbital uncertainty, thus making the dispersion due to the 
Yarkovsky the predominant source of uncertainty. 
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IAU Name GM 
No. (km 3 /s 2 ) 



1 


Ceres 


63.13 a 


4 


Vesta 


17.29 6 


2 


Pallas 


13.73 c 


10 


Hygiea 


5.78 a 


511 


Davida 


2.26 d 


704 


Inter amnia 


2.19 d 


15 


Eunomia 


2.10 d 


3 


Juno 


1.82 d 


16 


Psyche 


1.81 d 


52 


Europa 


1.59 d 


88 


Thisbe 


1.02 d 


87 


Sylvia 


0.99 d 


6 


Hebe 


0.93 d 


65 


Cybele 


0.91 d 


7 


Iris 


0.86 d 


29 


Amphitrite 


0.86 d 


532 


Herculina 


0.77 d 


324 


Bamberga 


0.69 d 


259 


Aletheia 


0.52 d 


11 


Parthenope 


0.39 d 


56 


Melete 


0.31 d 


14 


Irene 


0.19 d 


419 


Aurelia 


0.12 d 


63 


Ausonia 


0.10 d 


135 


Hertha25 


0.08 d 



2011 




Russell et al. 


(2012 





Konopliv et al 



d201lh . iCarrjl ((2012) 



Epoch 


2006 Oct. 23.0 


TDB 


Eccentricity 


0.1910634179(303) 




Perihelion distance 


0.7460469492(282) 


au 


Perihelion time 


2006 Jul. 7.8082161(101) 


TDB 


Longitude of Ascending Node 


204.4601395(387) 


deg 


Argument of perihelion 


126.3934237(392) 


deg 


Inclination 


3.33133146(361) 


deg 



Table 2: Apophis' orbital parameters. Numbers in parentheses indicate the 1-ct formal 
uncertainties of the corresponding digits in the parameter value. 



Solution 


N obs 


c 


c 






AC 


AC 






km 


km 


km 


km 


km 


km 


T 


433 


9487.57 


47683.39 


6.07 


285.74 


-1.24 


-24.15 


TR 


440 


9485.66 


47413.94 


5.95 


205.20 


0.67 


245.30 


TRO 


452 


9486.33 


47659.24 


5.75 


43.84 








All 


1603 


9482.98 


47684.20 


4.38 


34.24 


3.35 


-24.96 



Table 3: 2029 &-plane coordinates and uncertainties for different sets of observations: T is 
for only Tholen et al. observations, TR is for Tholen et al. plus radar observations, and 
TRO is for Tholen et al., radar, plus Pan-STARRS PS1 and Magdalena Ridge observations, 
ALL is the solution obtained by using all the available observations. A£ and A£ are with 
respect to the TRO solution. 
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Figure 2: Sampling of Apophis diameter D, Bond albedo A, density p, thermal inertia T, 
rotation period P ro t, and obliquity 7. 




Figure 3: Distribution of A2 and da/dt obtained from the assumed physical modeling. 
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Figure 4: 2036 resonant return, LOV, and keyhole on the 2029 6-plane. The impact cross 
section of the Earth is 2.16 i?®. 
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Figure 5: Keyhole map on 2029 6-plane. The minimum distance after 2029 is plotted as a 
function of £2029- The top panel shows the full uncertainty region (60000 samples), while 
the bottom panel is a zoom on the 2051 impactor complex (200000 samples). 
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Figure 6: Probability distribution on the 2029 6-plane due to the uncertainty in the orbital 
solution, the Yarkovsky effect, and both. The origin is set to the nominal C2029 = 47659.24 
km. Vertical bars correspond to keyholes and the height (in logarithmic scale) is propor- 
tional to the keyhole width. For clarity, only the keyholes with a width larger than 1 m 
are labeled with the year of impact. 
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Date width A( wrt TRO a IP 



TDB 




m 


km 








2042 Apr. 


13.5 


1.80 


-1557.3 


-4.80 


1.60 x 10- 


ii 


2041 Apr. 


15.3 


0.03 


-1549.4 


-4.78 


3.04 x 10- 


13 


2036 Apr. 


13.4 


616.16 


-1543.8 


-4.77 


7.07 x 10" 


-9 


2053 Apr. 


12.9 


3.09 


-1528.9 


-4.73 


4.91 x 10" 


11 


2044 Apr. 


13.3 


0.74 


-1523.1 


-4.71 


1.34 x 10" 


11 


2089 Oct. 


15.6 


0.05 


-788.7 


-2.38 


4.74 x 10' 


-9 


2060 Apr. 


12.6 


0.03 


-223.0 


0.08 


3.23 x 10- 


-8 


2065 Apr. 


11.8 


0.08 


-218.8 


0.10 


8.40 x 10- 


-8 


2078 Apr. 


13.8 


0.07 


-218.7 


0.10 


7.34 x 10- 


-8 


2091 Apr. 


13.4 


0.07 


-218.6 


0.10 


7.34 x 10- 


-8 


2103 Apr. 14.4 


0.18 


-212.5 


0.12 


1.82 x 10 


-7 


2077 Apr. 13.5 


0.06 


-212.3 


0.12 


6.05 x 10- 


-8 


2105 Apr. 


13.8 


0.03 


-212.3 


0.12 


3.03 x 10- 


-8 


2068 Apr. 12.6 


2.25 


-212.1 


0.12 


2.27 x 10 


-6 


2076 Apr. 13.0 


0.19 


-206.7 


0.13 


1.85 x 10 


-7 


2068 Oct. 


15.4 


0.14 


47.7 


0.54 


5.06 x 10- 


-8 


2105 Oct 


. 16.3 


0.38 


47.7 


0.54 


1.37 x 10 


-7 


2069 Oct 


. 16.0 


1.87 


671.6 


2.24 


2.39 x 10 


-7 


2069 Apr. 


13.1 


5.58 


1034.5 


3.53 


1.39 x 10- 


-8 


2077 Apr. 


13.6 


0.47 


1052.4 


3.59 


9.51 x lO- 


10 



Table 4: Risk assessment for Apophis. Columns are impact epoch, 2029 keyhole width 
and location, inverse of the normal cumulative density function, and impact probability. 
Bold lines correspond to a n impact probability larger than Sentry's generic completeness, 



i.e., a value of 8.6 x 10 8 ([Milani et al 



2005al) . 
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Assumption 



2036 IP 2036 AIP 2068 IP 2068 AIP 



Nominal 


7.07 x 10" 


-9 


Diameter 270 ± 60 m 


3.06 x 10" 


-6 


Unif. cos 7 


2.67 x 10- 


-9 


Symm. cos 7 


1.64 x 10- 


-9 


al = 0.6 


6.60 x 10- 


-7 


Vd = 200 


1.58 x 10- 


-9 


No = 400 


1.59 x 10- 


-8 



0% 


2.27 x 10~ 6 


0% 


+43169% 


1.68 x 10~ 6 


-26% 


-62% 


2.14 x 10~ 6 


-6% 


-77% 


1.06 x 10~ 6 


-53% 


+9242% 


2.38 x 10~ 6 


+5% 


-78% 


3.20 x 10~ 6 


+41% 


+125% 


2.02 x 10~ 6 


-11% 



Table 5: Variation of the impact probability for the 2036 and 2068 (AC = -212.14 km) 
keyholes for different assumption t he Yar kovsky effect modeling: nominal assumptions of 
SecG2 diameter from lDelbo et al.l (|2007al ); uniform distribution for cos 7; distribution for 
cos 7 symmetric to the one of Fig. [2] different variance of p; different means of do m the 
thermal inertia modeling. 




Figure 7: Arecibo delay-Doppler predictions for 2013-02-06.0 for different Yarkovsky sce- 
narios. The origin is set to (127.073088 s, -56954.4 Hz). 
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Scenario A 2 

1(T 15 au/d 2 



Date 
UTC 



Delay 



a) 2013-02-05.0 124.886779(1) 

a) 2013-02-06.0 127.073088(1) 

b) 25.0 2013-02-05.0 124.886797(1) 

b) 25.0 2013-02-06.0 127.073107(1) 

c) -25.0 2013-02-05.0 124.886760(1) 
c) -25.0 2013-02-06.0 127.073069(1) 



Table 6: Simulated February 2013 delay observations from Arecibo for different values of 
A2. The uncertainty is assumed to be 1 fis. 



Dates a A2 [10~ 15 au/d 2 ] 



Jan. 


2013- 


Jun. 2013 


47 


Sep. 


2013 - 


May 2015 


40 


Oct. 


2015 - 


Jan. 2016 


39 


Mar. 


2018 - 


May 2018 


32 


Dec. 


2018 - 


Aug. 2020 


16 


Nov. 


2020 - 


May 2021 


4 



Table 7: Uncertainty in the Yarkovsky effect determination as a function of increasing 
number of simulated optical observations weighted at 0.15" . For each interval we simulated 
one observation every two months. 
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